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Abstract 

We study a generalization of the fully overdamped Frenkel-Kontorova model in dimension n > 1. This model 
describes the evolution of the position of each atom in a crystal, and is mathematically given by an infinite 
system of coupled first order ODEs. We prove that for a suitable rcscaling of this model, the solution converges 
to the solution of a Peierls-Nabarro model, which is a coupled system of two PDEs (typically an elliptic PDE 
in a domain with an evolution PDE on the boundary of the domain). This passage from the discrete model to 
a continuous model is done in the framework of viscosity solutions. 
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1 Introduction 

In this paper we are interested in two models describing the evolution of defects in crystals, called 
dislocations. These two models are the Frenkel-Kontorova model and the Peierls-Nabarro model. The 
Frenkel-Kontorova model is a discrete model which describes the evolution of the position of atoms in 
a crystal. On the contrary, the Peierls-Nabarro model is a continuous model where the dislocation is 
seen as a phase transition. The main goal of the paper is to show rigorously how the Peierls-Nabarro 
model can be obtained as a limit of the Frenkel-Kontorova model after a suitable rescaling. 

1.1 Peierls-Nabarro model 

Let us start to present the Peierls-Nabarro model. We set 

n = {x = (x 1 ,...,x n ) E R n , x n > 0}, 
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and for a time < T < +00, we look for solutions u° of the following system with /3 > 0: 

Pu° t {x,t) = Au°(x,t) (x, t) e O x (0, T) 

<9n° (1-1) 
u° t (x, t) = F(u°{x, t)) + -k—{x, t) (x, t) G x (0, T), 

Cx n 

where the boundary 9f2 is defined by: 

<9f} = {x = (xi,...,x n ) G M n , x n = 0}, 

and the unknown u°(x,t) G R is a scalar- valued function with the initial data 

f x G H := <9ft U Q if /3 > 
u°(x,0) =u (x) for < 1.2 

v ' ' v ; 1 x G dft if p = 0. v ; 

We assume the following conditions on the function F : R — > R and on the initial data ito : 

f n G ^ 3 '°°(0) if p > 
FG^ 2 '°°(M) and ^ , 7 ^ 1.3) 

\ n G VF 3 ' 00 ^) if /3 = 0. 

For simplicity we have taken a high regularity on the initial data, but this condition can be weaken. 
The classical Peierls-Nabarro model corresponds to the case (5 = with a 1-periodic function F (see 
Section 2). For a derivation and study of the Peierls-Nabarro model, we refer the reader to [22] (and 
[30, 28] for the original papers and [29] for a recent review by Nabarro on the Peierls-Nabarro model). 
Let us mention that a physical and numerical study of the evolution problem (1.1) has been treated 
in [27]. 

In this paper we consider the general case (5 > 0, since the case (3 = is natural in the Frenkel- 
Kontorova model (see Section 2), and mathematically the case P > is easier to study. In the special 
case P = 0, this problem can be reformulated (at least at a formal level) as a nonlocal evolution 
equation written on d£l in any dimension. See in particular [18] for such a reformulation in dimension 
n = 2 and a limit of the Peierls-Nabarro model to the discrete dislocation dynamics (see also [26] for 
a homogenization result of the Peierls-Nabarro model). We refer the reader to [10] for the study of 
stationary solutions of our model when F = —W'. 

1.2 Prenkel-Kontorova model 

We now present the Frenkel-Kontorova model. This is a discrete model which contains a small scale 
e > that can be seen as the order of the distance between atoms. We set the discrete analogue of 
O: 

fl £ = (£-L) n - x x e(N \ {0}) 

with discrete boundary 

dn £ = (eZ)™" 1 x {0}, 



and we set 



n £ = n e udn E . 



We look for solutions u £ (x, t) with x G Q £ , t > 0, of the following system for P > 0: 

f Pu £ t (x,t) = A £ [u £ ](x,t), (x,t) en £ x (0,T) 

I u £ (x,t) = F(u £ (x,t)) + D £ [u £ ](x,t), (x,t) £dn £ x (0,T), 



(1.4) 
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with initial data 



u £ (x, 0) = uq(x) for 



if ft > 



xedn £ if ft = o, 



(1.5) 



where we have used the following notation for the discrete analogue of A and -Jj— respectively: 

A £ [u £ ]{x,t) = \ (u £ (x + ey,t) -u £ {x,t)) 



y&Z n ,\y\=l 



D £ [u £ ](x,t) 



(u £ (x + ey,t) -u £ (x,t)). 



;i.6) 



y&z n ,\y\=i,yn>o 



In this model, the function u £ (x, t) describes the scalar analogue of the position of the atom of index x 
in the crystal. The classical fully overdamped Frenkel-Kontorova model is a one dimensional model. 
It corresponds to the particular case n = 2 where u £ (x,t) = if x E £1 £ , and therefore only the 
u £ (ex,t) for x E Z x {0} can be non trivial (we refer the reader to [9]). Remark that in the case 
ft > 0, we can consider classical solutions using the Cauchy-Lipschitz theory, while at least in the 
case ft = 0, it is convenient to deal with viscosity solutions (see Section 3). 

More generally, Frenkel-Kontorova are also used for the description of vacancy defects at equi- 
libruim, see [20]. See also [19], where the authors study the problem involving a dislocation inside 
the interphase between two identical lattices. Their model corresponds to our model (1.4) at the 
equilibrium with F = —W where the potential W is a cosine function. For other 2D FK models, see 
[12, 13]. For homogenization results of some FK models, we refer the reader to [17]. 



1.3 Main result 

Our main result is the following theorem which establishes rigorously (for the first time up to our 
knowledge) the link between the two famous physical models: the Frenkel-Kontorova model and the 
Peierls-Nabarro model. 



Theorem 1.1 (Existence, uniqueness and convergence) 

Let e > 0, ft > and < T < +oo. Under the condition (1.3) there exists a unique discrete viscosity 
solution u £ of (1.4)-(1.5). Moreover, as e — > 0, the sequence u £ converges to the unique bounded 
viscosity solution u° of (1.1)-(1.2). The convergence u £ — > u° has to be understood in the following 
sense: for any compact set K C O x [0, T), we have: 

W uS ~ u llL°°(ii:n(o e x[o,T))) ~^ as e ~^ °- 

A formal version of this result has been announced in [16]. We refer the reader to Section 3 for 
the definition of viscosity solutions. The proof of convergence is done in the framework of viscosity 
solutions, using the half-relaxed limits. The uniqueness of the limit u° follows from a comparison 
principle that we prove for (1.1), using in particular a special test function introduced by Barles in 
[5]. Most of the difficulties arise here from the unboundedness of the domain, and in the non standard 
case ft = 0, where the evolution equation is on the boundary of the domain instead of being in the 
interior of the domain. Let us notice that because system (1.4) can be seen as a discretization scheme 
of system (1.1), then Theorem 1.1 can be interpreted as a convergence result for such a scheme. In 
the literature on numerical analysis of finite difference schemes, we find other results of convergence 
for different equations (see for instance [8] and [7]). Let us mention that an estimate on the rate of 
convergence of u £ to u° is still an open question for our system. 
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Remark 1.2 From the proof of Theorem 1.1, it is easy to see that the result still holds true if F 
is replaced by a sequence (F e ) £ of functions which converges in W 2 '°°(R). For such an example of 
application, see (2.3) at the end of Section 2; see also [16]. 



1.4 Organization of the paper 

This paper is organized as follows. In Section 2, we present the physical motivation to our problem. 
In Section 3, we present the definitions of viscosity solutions for the discrete and continuous problem 
(1.4)-(1.5) and (1.1)-(1.2) respectively. Section 4 is dedicated to construct uniform barriers of the 
solution u £ of (1.4)-(1.5). Using those barriers, we prove in Section 5 the existence of a solution for 
the discrete problem (1.4)-(1.5). Section 6 is devoted to the proof of Theorem 1.1. Sections 7 and 8 
are respectively devoted to the proofs of the comparison principle for the discrete problem (1.4)-(1.5) 
and the continuous problem (1.1)-(1.2), that were presented in Section 3 and used in the whole paper. 
Finally, in the Appendix, we give a convenient corollary of Ishii's lemma, that we use in Section 8. 



2 Physical motivation 

2.1 Geometrical description 

We start to call (e\, &2% e z) an orthonormal basis of the three dimensional space. In the corresponding 
coordinates, we consider a three-dimensional crystal where each atom is initially at the position of a 
node I of the lattice 

A = Z 2 x A + 

This lattice is simply obtained by a translation along the vector ^3 of the lattice Z 3 , and will be 
more convenient for the derivation of our model. We assume that each atom / 6 A of the crystal has 
the freedom to move to another position I + u\ei where ui E R is a unidimensional displacement. 
In particular, we will be able to describe dislocations only with Burgers vectors that are multiples 
of the vector e 2 (see [22] for an introduction to dislocations and a definition of the Burgers vector). 
Moreover we introduce the general notation 

7=(/ / ,-/ 3 ) for all /= (I',I 3 ) eZ 2 x (- + Z) with /' = (I u I 2 ) € Z 2 
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and will only consider antisymmetric displacements uj, i.e. satisfying 

uj = —ui for all I G A. (2.1) 
We also introduce the discrete gradient 

(vV)j 

Generally, the core of a dislocation is localized where the discrete gradient is not small. In our model, 
screw and edge dislocations will be represented as follows. 



Ul+ei 


- ui 


Ul+e 2 


- ui 


u I+e 3 


- Ui 
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Screw dislocation 

For / = (Ii, I2, 13), we can consider a dislocation line parallel to the vector e2, with a displacement 
ui independent on I2, i.e. such that uj = s(Ii,Is). We will assume for instance that 

(s(i"i,+i) = -s(Ii,-~) ->• as Ii ->■ -00 
s(Ji, +i) = -s(Ji, -5) -)• 5 as Jx-^+oo. 
Moreover, if there is no applied stress on the crystal, then it is reasonable to assume that 

dist ((V d it)j,Z 3 ) -> as | (J l9 1 3 )| -> +00 

which means that the crystal is perfect far from the core of the dislocation. 
Edge dislocation 

For I = (Ii, I 2 , 13), we can consider a dislocation line parallel to the vector ex, with a displacement 
uj independent on Ix, i.e. uj = e{l2,Is)- We will assume for instance that 

e(I 2) +|) = — e(J 2 , — \) — > as I 2 ->■ — 00 

e(J 2 , +5) = — e(I 2 ,— |) ->■ i as J 2 -> +00. 
Moreover, if there is no applied stress on the crystal, then it is reasonable to assume that 

dist ((V^) 7 ,Z 3 ) ->-0 as \(h,h)\ -> +00. 

Remark that this model can also describe more generally curved dislocations, which are neither screw, 
nor edge, but are mixed dislocations. 



2.2 Energy of the crystal 

We assume that each atom / is related to its nearest neighbors J by a nonlinear spring, whose force is 
derived from a smooth potential Wjj. Then formally the full energy of the crystal after a displacement 
u = (uj)i 6 a is 

E{u) = \ W L j{ Ul - uj). 

/,JeA, |/-J|=i 



We assume that the dislocation cores are only included in the double plane 1 3 = We also assume 
that 

/ eW(a) if |7 3 | = I J 3 | = \ 
Wlj{a) = \w(a) if |/ 3 |/I or Val^i 

Here e > is a small parameter. This means that the springs lying between the two planes I3 = ±^ 
are very weak in comparison to the other springs. This will allow the core of the dislocation to spread 
out on the lattice as e — > and will allow us to recover the Peierls-Nabarro model in this limit, after 
a suitable rescaling. 

Notice that, to be compatible with the antisymmetry (2.1) of the crystal, it is reasonable to 
assume that 

W(-a) = W(a) for all a £ M. 

Recall that each line of atoms 1° + Ze2 only contains atoms of identical nature (i.e. surrounded by 
a similar configuration of springs). Therefore for ui = s(Ii,Is), a lattice of atoms at the position 
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I + Ufe2 and another lattice of atoms at the position / + u\&2 + k(ii,i 3 ) e 2 with arbitrary fe^ j 3 ) £ Z 
are completely equivalent. Therefore their energy should be the same, and it is then natural to 
assume that the potential W is 1-periodic, i.e. satisfies 

W(a + 1) = W(o) for all a G R. 

We will also assume that non-deformed lattices minimize globally their energy, i.e. satisfy 

W(Z) = and W > on R\Z. 

Let us now assume that a constant global shear stress is applied on the crystal such that it creates 
a global shear of the crystal in the direction e<i with respect to the coordinate I3. This means that 
there exists a constant r £ R such that 

dist ((V d u)/ - re 3 ,Z 3 ) -»• 

as / goes far away from the core of the dislocation. Remark that this assumption is compatible with 
the antisymmetry (2.1) of u. If there is no dislocations, we can simply take for instance uj = r/3. 

2.3 Fully overdamped dynamics of the crystal 

The natural dynamics should be given by Newton's law satisfied by each atom. This dynamics is very 
rich and for certain shear stress r, it is known (see [14, 21]) in 2D lattices that certain edge dislocations 
can propagate with constant mean velocity. This is due to the fact that part of the energy is lost by 
radiation of sound waves in the crystal. This phenomenon is similar to the effective drag force created 
by the surrounding fluid on a boat or an air plane. The resulting behaviour is a kind of dissipative 
dynamics that we modelise here by a fully overdamped dynamics of the crystal. See [23, 25, 24]) 
for a fundamental justification of overdamped type dynamics based on explicit computations in a 
ID Hamiltonian model. For general physical justifications of the dissipative effects in the motion of 
dislocations, see also [22, 1]. 

We recall that we assume that the dislocation cores are only contained in the double plane I3 = =t|. 
For this reason, we will artificially distinguish the dynamics inside this double plane and outside this 
double plane. We consider the following fully overdamped dynamics (where the velocity of each atom 
is proportional to the force deriving from the energy) which is written formally as: 

aiiii = —V Ul E(u) with OLi = i ^ ["j 3 ! , f (2-2) 

for some constant a > 0. Here iii denotes the time derivative of the displacement uj and 

V Ul E(u):= Yl W'jAm-uj). 
JeA, \J-i\=i 

Remark in particular that this dynamics preserves the antisymmetry (2.1) of u. For slow motion of 
dislocations (i.e. small velocity with respect to the velocity of sound in the crystal), it is reasonable 
to assume that a = 0, i.e. the lattice is instantaneously at the equilibrium outside the double plane 
h = ± l 2 - 

We also assume that the potential is harmonic close to its minima, i.e. satisfies 

a 2 1 
W(a) = — for \a\ < 5 < - 
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and assume that the strain is small enough in the crystal outside the double plane 1% = ±^> i- e - we 
assume that 

\uj — uj\ < 5 for any J, J € A such that | J — I\ = 1 and | _?3 1 7^ — or | J3I 7^ — . 

This assumption allows us, in the region I3 ^ to consider forces that can be expressed linearly 
in terms of the displacement. We also set 

t = ea for some a G R and vj = m — eal^. 

Using the antisymmetry of the solution (2.1), we deduce that we can rewrite the dynamics (2.2) with 
a = as 

= ^2 ( VJ ~ v ri for /3 > 2 

JeA, |J-/|=i 

< 

vi = ea — eW'(2vi + ea) + {vj — vi) for ^3 = -. 

,/eA, |J-/|=1, J 3 >| 

Then we see that 

U £ (x, t) = V^x + l-j 

solves (1.4) with [3 = and F replaced by 

F £ (a) = a- W'{2a + ea). (2.3) 

Remark 2.1 At the level of modeling, we could also consider more general lattices than Z n with 
more general nearest neighbors interactions, but this case is not covered by the result of Theorem 1.1 
and would require a specific work. 




3 Viscosity solutions 

In this section we present the notion of viscosity solutions and some of their properties for the 
discrete problem (1.4)-(1.5) and then for the continuous problem (1.1)-(1.2). For the classical notion 
of viscosity solutions, we refer the reader to [3], [4], [15]. 

3.1 Viscosity solutions for the discrete problem 

Before stating the definition of viscosity solutions for the discrete problem (1.4)-(1.5), we start by 
defining some terminology. Let e,5 > and < T < +00. Given a point Pq = (xo,to) G Q x [0, T), 
the set yV|(t ) cffx [0,T) is defined by: 

A/|(P ) = {P = (V,t) en £ x [0,T);\y - x \ < e and \t-t \<5}. (3.1) 

The spaces USC(Q x [0,T)) and LSC(Q x [0,T)) are defined respectively by: 

USC(Tl £ x [0,T)) = {u defined on Q £ x [0,T); Vx G ff , u(x, • ) G USC{[0,T))} (3.2) 
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and 

LSCiSf x [0,T)) = {u defined onf x [0,T); Vxeff, u(x, • ) G LSC([0,T))}, (3.3) 

where USC([0,T)) (resp. LSC([0,T))) is the set of locally bounded upper (resp. lower) semi- 
continuous functions on [0, T). In a similar manner, we define the space C k (Q £ x [0,T)), k G N, 

C fc (lf x [0,T)) = {u defined on if x [0,T); Vx G Q £ , u(x, • ) G C fc ([0,T))}. (3.4) 
Next, given f3 > 0, we present the following definition: 

Definition 3.1 (Viscosity sub/super solutions) 

1. Viscosity sub-solutions. A function u G USC(il x [0, T)) is a viscosity sub-solution of problem 
(1.4)-(1.5) provided that: 

' x G IF if > 
xedQ £ if j3 = 0, 

(ii) for any ip G C 1 (r2 x [0, T)) ; ifu — ip has a zero local maximum at a point Pq = (xq, to) € Q x [0, T) 
(i.e. 3 5 > such that VP G A/|(P ), we /wue (u - y?)(P) < (n - <p)(Po) = 0) then 



(i) n(x, 0) < uq(x) for 



for P G il £ x (0, T) if P>0 
for P G ft e x [0, T) if P = (3.5) 
^(P ) < P(v?(P )) + Z^MtPo) /or P G 0fi e x (0,T). 



M(P )< A e M(P ) 



2. Viscosity super-solutions. A function u G LSC(il x [0,T)) is a viscosity super- solution of 
problem (1.4)-(1.5) provided that: 

' x G Tf if j3 > 
i£ffi £ i/ /? = 0, 

(ii) for any ip G C l (il x[0,T)),ifu — ip has a zero local minimum at a point Pq = (xQ,to) G il x [0, T) 
(i.e. 3 5 > such that VP G A/|(P ), we We (n - ip)(P) > (u £ - <p)(P ) = 0) then 



(i) u(x,0) > «o(x) /or 



/or P G ft e x (0, T) if > 
/or P G ft e x [0, T) i/ = (3.6) 
I ^(^o) > P(^(Po)) + ^ £ M(^o) for P G dO £ x (0,T). 



/^(P ) > A £ M(^o 



3. Viscosity solutions. A function u G C°(il £ x [0, T)) is a viscosity solution of problem (1.4)-(1.5) 
if it is a viscosity sub- and super-solution o/ (1.4) -(1.5). 

Definition 3.2 VFe say i/iai u is a viscosity sub-solution (resp. super- solution) of (1.4) if u only 
satisfies (ii) in the Definition 3.1. 

Remark 3.3 Note that a function u is a viscosity subsolution of (1.4) if it satisfies the viscosity 
inequalities on if x (0,T) if p > and on (if x (0,T)) U (il £ x {0}) i//3 = 0. 

Theorem 3.4 (Comparison principle in the discrete case) 

Letu G U SC(if x [0, T)) (resp. v G LSC(ifx [0, T))) 6e a viscosity sub-solution (resp. super- solution) 
for the problem (1.4)-(1.5) where < T < oo, suc/i i/iai ||^|lioo(Q e x [ y)), IMlL°°(n £ x[o T)) < Then 

u(x,t) <v(x,t) for all (x, i) £(f x [0, T). 
Theorem 3.4 will be proved in Section 7. 
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3.2 Viscosity solutions for the continuous problem 

Similarly as in subsection 3.1 we denote by USC(Q x [0,T)) (resp. LSC(Q x [0, T))) the set of 
locally bounded upper (resp. lower) semi-continuous functions on x [0, T). We also denote by 
C k (Q x [0, T)) the classical set of continuously /c-differentiable functions on 17 x [0, T). 

Definition 3.5 (Viscosity sub/super solutions) 

1. Viscosity sub-solutions. A function u G USC(Q x [0, T)) is a viscosity sub-solution of problem 
(1.1)-(1.2) provided that: 

' H x {0} if p > 
dU x {0} if (3 = 0. 

(ii) /or any <p £ C 2 (r2 x [0, T)), ifu — ip has a zero local maximum at a point Pq = (xq, to) G ^ x [0, T) 
t/jen 

/or P G x (0, T) if (3 > 



(i) u <uq on 



p<p t (P ) < A^(Po) 



/orPoG^x [0,T) i/ = 



min|M(P ) - A^(P ) , Pt(fb) - P(^(Po)) - J^( P o) j < /or P € fl« x (0,T). 

(3-7) 

2. Viscosity super-solutions. A function u £ LSC(Vt x [0, T)) is a viscosity super- solution of 
problem (1.1)-(1.2) provided that: 

' U x {0} if P>0 
dQ x {0} »/ /3 = 0. 

(ii) /or any <p> G C 2 (f2 x [0, T)) ; ifu — p has a zero local minimum at a point Pq = (xq, to) £ x [0, T) 
t/ien 

/or P e H x (0, T) if > 



(i) u> uq on 



P<p t (P ) > A<p(P ) 



forPoeflx [Q,T) if P = 



max j/%(P ) - A^(P ) , tp t (P ) - F(^(P )) - ^r^o) j > for P e dfl x (0,T). 

(3-8) 

3. Viscosity solutions. A function u £ C°(Q x [0, T)) zs a viscosity solution of problem (1.1)-(1.2) 
if it is a viscosity sub- and super- solution of (1.1)-(1.2). 

Theorem 3.6 (Comparison principle in the continuous case) 

Letu € USC(£lx [0, T)) (resp. v £ LSC(£l x [0,T))) be a viscosity sub-solution (resp. super- solution) 
of problem (1.1)-(1.2) where < T < oo, suc/i that IM|£oom x ro t))> ll u llx,°=>(T5x[o t)) ^ +°°- 

u < w on H x [0, T) 

This theorem will be proved in section 8. 

Remark 3.7 Remark that in the case (3 = 0, the problem can be reformulated (at least for smooth 
solutions) as a nonlocal evolution equation written on the boundary d£l. See for instance the work 
[11] on the relation bewteen fractional Laplacian and harmonic extensions. Note that there is also a 
viscosity theory for nonlocal operators (see [6]). 

4 Construction of barriers 

This section is devoted to the construction of barriers for the solution u £ of (1.4)-(1.5) for all P > 0. 
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4.1 Discrete harmonic extension 

Here we construct the discrete harmonic extension of uq on f2 e , i.e. we prove the existence of a 
solution of the following problem 

f - A 1 ^] =0, inO 1 

I Uq = Uq, On Oil , 

where A 1 = A e for e = 1 is defined in (1.6). We note here that, for the sake of simplicity, we have 
taken e = 1, while otherwise it can be treated in the same way (or simply deduced by rescaling). We 
will say that a function u is a sub-solution (resp. super-solution) of (4.1) if the symbol "="in (4.1) is 
replaced by the symbol "<"(resp. ">"). 

Lemma 4.1 If uq G L^^dQ 1 ), then there exists a unique (viscosity) solution o/(4.1) on Q, . 
Proof. 

Step 1: Existence. 

Define 

u = sup no and u = inf uq , 

an 1 dn 1 

then u (resp. u) is a sub- (resp. super-) solution of (4.1). Consider now the set 

S = |u : fJ 1 — > R, subsolution of (4.1) s.t. u<u< n| , 
and in order to use Perron's method, we define 

Uq := sup?/. 

ues 

Using the fact that the maximum of two sub-solutions is a sub-solution, it is possible to show that Uq 
is a sub-solution of (4.1). Hence, it remains to prove that Uq is a super-solution for (4.1). Suppose 
that there exists xq € Q 1 such that 

- A^K^o) < 0, ifxoG^ 1 



or 



k Uq(x )<u (x ), iixoGdn 1 . 
Then we construct a sub-solution w G S such that w > Uq at xq. Two cases are considered. 
Case -A 1 [u D ](x ) < for x G ft 1 . Let w be defined by 



W ( X ) = { ..Df_ N , 1 / 



Uq (x) if X ^ Xo 

Uo(x ) - 

We check that u<w<u, and we compute 



u o ( x o) + 7T A i u o }( x o) if x = x Q . 
In 



>A 1 [«5 , ](x)>0 if x^x 
if x = xq, 
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where we have used the fact that w > Uq to deduce the inequality. Therefore, we have w E S and 
w(xq) > Uq(xq) which is a contradiction. 

Case Uq ) (xo) < uo(xo) for xo E dfl 1 . Let w be defined by 

, , / Uq(x) if x ^ x , 
[ uq(xq) if re = x . 

By construction, we have u < w <u. Moreover, as w > Uq it follows that A 1 ^] > A^Iuq] > on 
Q 1 . Therefore w£5 and w(xq) = uq(xq) > Uq(xq); this implies a contradiction. 

Step 2: Uniqueness. 

We simply adapt case 1 i) of the proof of the comparison principle (Theorem 3.4), given in Section 
7. □ 

4.2 Continuous harmonic extension 

For the case /3 = 0, we need to consider the harmonic "extension" Uq : — > M of the initial data uo- 
This function is the solution of 

" -Atig = on n 

Uq = Uq On Oil. 

This section is devoted to recall the existence of the continuous harmonic extension and to show 
some of its properties. Let z = (V, z n ) E Q where z' is identified to an element of <9f2 and z n > 0, 
and let H be a new function defined by 

where ui n is the measure of the (n — l)-dimensional sphere in M n . Now, we define Uq by 

u c q{x) := / H(x' - z', x n )uo(z') dz' for all x = (x' , x n ) E W 1 ' 1 x (0, +oo). (4.4) 
Jan 

Note that H(x' — z',x n ) is the Poisson kernel. We have the following three lemmas: 
Lemma 4.2 (Existence of a continuous harmonic extension) 

If uq(z') is bounded and continuous for z' E d£l, then the function Uq(x) defined by (4.4) belongs to 
C°°(f2) and is harmonic in Q and extends continuously to O such that Uq = uq on d£l. 

Proof. See [2, Section 7.3, p. 129]. □ 

Lemma 4.3 (Estimate on the harmonic extension) 

IfuQ E W 3 '°°(dn), then we have 

C C 
\Du C q(x)\ < , \D 2 u C q(x)\ < j for allx = (x',x n ) E R^ 1 x (0,+oo). (4.5) 

1 ~\~ Xn 1 ~\~ X n 

Proof. Since no is bounded, then by adding an appropriate constant to uq we can always assume, 
without loss of generality, that uq > 0. From (4.4) and from the following property (see [2, p. 126]) 

/ H(x — z'i x n ) dz = 1 for all x = (x \ x n ) E 17, 
Jdn 
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we have: 



< / H(x - z', x n )\uo(z')\ dz' < \\uq\\l°° / H(x - z ,x n ) dz 
'an Jan 



and then 

IKIU°° < |holk°°- (4.6) 
Next, for x = (x',x n ) £ W 1 ' 1 x (0, +00), using (4.3) and (4.4), we have, on the one hand with \x = — : 



1 du c n 



u Q (z') 



Hdx n J dn (x 2 + {x> - z') 2 ) n / 2 



dz' 



nx 2 l uo(z') 
n Wn + - z') 2 ) n/2+1 



dz' , 



(4.7) 



then we obtain 



du c n 



dx r , 



-u,il°° w jMk°° < _c 



where we have used (4.6). On the other hand, we have: 

x n u {z')(x' - z') 



1 8uq 
fi dx' 

then, using Young's inequality: 



-n 



dQ (x* + (x> - z'YY^ 



x n \x' - z'\ < ^-(x z n + {x' - z' f) 



dz' 



(4.8) 



(4.9) 



we conclude that 



Therefore 



duf 



dx' 



< 



Uq(z') 



dz'< 



n 

n {xl + {x> - z'YYl* "~ 2x~ r , 

c 



< 



c 



\Du c (x)\ < 

In order to prove the second inequality in (4.5), we differentiate (4.7) and (4.8), we get 



1 d 2 u^ f u (z')[-nx n I'} 
H dx'dx' ~ J m (x 2 n + (x' - Z ')2)n/2+l 



dz' + 



x n u (z')(n)(n + 2)(x' - z') ® jx' - z') , 
an [xi 



[Xt + (x' - Z 'V)n/2+2 



where I 1 is the identity matrix of W 1 , then, as l/(cc„ + (x' — z') ) < we obtain 



dh 







dx'dx' 



< n 



u (z')x n 



av. 



[xl + {x> - z') 2 y/ 2 xl 



dz'+n(n + 2) 



uo(z')x n (x' - z') 2 



9n (xl + (x' -z'yy/Hxl + (x' -z'Y)' 



dz'. 



Thus, using the following inequality 



we deduce that 



(x'-z'f 1 



(x 2 +(x'-z') 2 ) 2 - xV 



d 2 u? 







dx'dx' 



C 



(4.10) 
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Moreover, we have 



1 d 2 u c 
p dx'dx r . 



u (z')(-n)(x' - z') 



dz' + 



' dn (* 2 + (x> - z') 2 )™/ 2+1 
J_ f u (z')x n (-n)x n (x' - z') 

4 J d n (xl + (x> - z<yyi\xi + [ x > - Z 'Y) 

Then, using (4.6) and (4.9), we infer that 



u (z')x 2 n (n)(n + 2)(x' - z') , 

Q {Xl + (X> - Z 'f)n/2+2 Z > 



dz' + 



u (z')xj(n)(n + 2)(x' - z') , 
(xl + (x' - z') 2 ) n / 2 + 2 Z ' 



d 2 u c 



dx'dx r 



- 9_ 

x 2 * 



Next since, 



1 d 2 U C n 



-3n 



fan (xl 



u (z')x r , 



p dx n dx n -"Jon (xl + (x> - z') 2 )"/ 2+1 
we use similar arguments as above in order to obtain 



dz' + n(n + 2) 

Jan 



u (z')xl 



(x 2 n + (x> - Z ') 2 y/ 2 + 2 



dz', 



d 2 u c 



dx n dx n 



- 9_ 

xQ 



Notice that this last inequality also follows from (4.10) joint to the harmonicity of Uq. 
Finally, we use Schauder's estimate near the boundary d£l : 

|""olc 2 ^(Si(a:',o)no) - ^ (\^ u o\c a (B 2 (x' ,o)nQ) + \ u o\c 2 - a {B 2 (x> ,o)ndn)) = C| n olc* 2 ^(B 2 (a:',o)n9n) 5 

for all x' G IR™" 1 , where C a and C 2 ' a , a E (0,1), are the Holder spaces of order a and 2 + a, 
respectively. The ball B r (x',0), r = 1,2, stands for the ball of center (x',0) and radius r. We finally 
conclude that: 

sup \uo\ C 2,<*(Bi(x' o)nn) — ^ SU P \ u o\c 2 ' a (B 2 (x',o)ndn): 



and the result follows. 

Lemma 4.4 (e-uniform bound) 

Under the assumption (1.3) we have 



□ 



\D £ K](x)\ < 



Ci 



|A e [^](*)| < 



for all x = (x', x n ) G Sl £ 



(4.11) 



1 ~\~ x n 1 ~\~ xl 

where C\ > is a positive constant independent on e. 

Proof. Let us set p~ = and pf = 1 otherwise. Then using Taylor's expansion, we have for x G Q £ : 

n n „i 

^ E K]W:=EE^*±^)-«§(i)) =EE^ / dtDul(x±ete l )(±ee l ) 1 
± i=l ± i=i ^ 

where e$ = (0, . . . , 0, 1, 0, . . . , 0) is the unit vector in M n with respect to the i component. Then, using 
(4.5), we get 

\D £ [u c ](x)\ < (2n- 1) sup / dtpf\Du c (x±etei)\ < - — — . 

Ki<n Jo l + Xn 
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Then, in order to terminate the proof, it is sufficient to write A £ [uq] using Taylor's expansion as 

n 

e 2 A^](x) := ^^(«S(x± eei )-«S(x)) 



± i=i 



^2^2^±ee l Vu c (x) + dt dsD 2 Uq(x ± esej)(±eej)(±eej)^ (4.12) 



e 2 V]Yj / di / dsZ? 2 Uo(x ± esej) • (ej, e,). 
± i =i- / o ./o 



Thus, thanks to (4.5), we finally obtain: 



|A £ [no](x)| < 2n sup / dt [ ds\D 2 ul(x ± ese;)| < - 



+ x 2 



/or (x,i) 6 fi e x [0,T), (4.13) 



and the proof is done. □ 
4.3 Uniform barriers for {3 > 

In this subsection, we show uniform barriers for the solution u £ of (1.4)-(1.5) for all (3 > in the 
special case where uq = Uq. 

Proposition 4.5 (Uniform barriers in e for all (3 > for uq = Uq) 

Under assumption (1.3), there exists a constant C > independent on e > 0, f3 > and T sitc/i i/iai 
if 

H+(x, t) := u§(x) + C( yi + x» - 1 + t) 
H~(x, i) := -Uq(x) - C(yjl + x„ - 1 - t) 

thenu + (resp. vT ) is a supersolution (resp. subsolution) of (1.4)-(1.5). However, if u £ is a bounded 
viscosity solution of (1.4)-(1.5), i/ien u>e /lave 

IT < u £ < H + on ff x [0, T), (4.14) 

and moreover: 

\u e (x,t)\ < \\u \\l°° /or (x, t) 6 Q £ x [0, T). 

Proof. 

Step 1: sub/supersolution property 

We first check that u + is a super-solution of (1.4)-(1.5). Indeed, 

• If (x,t) G O, 6 x {0}, then H + (x,i) > Uq(x) = Uo(^)- 

• If (x,t) £ dn e x (0,T), then we take C > such that 

C > F(H+(x,t)) +D e [u[j](x) + C J D e [Vl + x„] H 4 + (x,t) > F(u + (x,t)) + D £ [u + ](x,t). 

Using the fact that D £ [^l + x n J(x n = 0) < |, we see that such a constant C exists. 

• If (x,t) £ £l £ x [0,T), then we can use Lemma 4.4. We get 



4(1 + x„) 3 / 2 ' 
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where the last inequality is true for C > large enough. Moreover, by repeating the same computa- 
tions as in (4.12), we easily get with f(a) = yl + a: 

A £ [VT+^} = f dt f dsY^f"{x n ±es) < f dt f ds.2.f"(x n ) < f"(x n ) = - 1 

Jo Jo ± Jo Jo 4(1 + x n yi 2 

(4.15) 

where we have used the concavity of the function /". Then we have 

C 



4(1 + / 

hence 

puf(x,t) = PC > > A £ [u + }(x,t). 

As a consequence, we finally deduce that u + is a super-solution. In a similar way, by taking in 
addition the following assumption on C : 



C > -F(u-(x,t)) - D £ [uq](x) + CD e y\ + x n ] 
we can prove that u~ is a sub-solution. 
Step 2: bounds on u e 

Let us call v(x,t) := H^oHl/ 00 + ^II-^IIl 00 - It is easy to check that v is a super-solution of (1.4)-(1.5). 
Then 

min(u + , v) 

is still a super-solution and is bounded. Therefore, if u £ is a bounded viscosity solution of (1.4)-(1.5), 
we can then apply the comparison principle (Theorem 3.4) to conclude that 

u £ < min(u + , v). 

This shows the upper bounds. For the lower bounds, we proceed similarly with max(u~, —v). □ 

4.4 Barriers for (3 > 

We have the following 

Proposition 4.6 (Barriers for j3 > 0) 

Under the assumptions of Theorem 1.1, for every (3 > 0, there exists a constant Cp > such that for 
all e, if u £ is a bounded viscosity solution of (1.4)-(1.5), then we have 

\u £ (x,t) -u (x)\ <tCp for (x, i) 6 ff x [0, T). 

Proof. 

Using in particular the fact that F is bounded and that uq G W 2 '°°(Q), we simply check that 
uo(x) + tCp is a super-solution for a suitable constant Cp> and apply the comparison principle. 
We proceed similarly for sub-solutions u$(x) — tCp. □ 



15 



5 Existence and uniqueness of a solution for the discrete problem 

The aim of this section is to prove the existence of solutions of problem (1.4)-(1.5). Cauchy-Lipschitz 
method is the main tool used to prove the existence of solutions for f3 > 0, while in the case (3 = 
we need barriers to prove the existence. 

Theorem 5.1 [Existence and uniqueness, (3 > 0) 

IJuq and F satisfy (1.3) and f3 > 0, then there exists a unique bounded solution u^' £ G C^{Vf x [0,T)) 
of problem (1.4)-(1.5). 

Proof. The proof is done using the classical Cauchy-Lipschitz theorem. Let B := L°°(f2 e ) be the 
Banach space with the norm 

\\u\\b ■= sup I , for every u G B, 
xeTi e 

and J- : B — > B be the map defined, for every u G B and x G f2 e , by 



F(u(x)) + D £ [u](x), x edn 

where A £ [u](x) and D £ [u](x) are defined as in (1.6), dropping the variable t on both sides of (1.6). 
Then, for every u,v G B and x 6 f2 , we have two cases; either x G D, £ , and hence we obtain 

\\F[u](x) - F[v](x)\\ B < \\±A £ [u-v](x)\\ B 

An 

or x G <9f2 £ , then: 



||J"M(x) -T\v\(x)\\ B < \\D £ [u-v}(x) + \F(u(x))- F(v(x))\\ 

[2{2n- 1) \ „ 



In all cases we conclude that J- is globally Lipschitz continuous. By the Cauchy-Lipschitz theorem, 
we get the existence and uniqueness of a solution u 13,6 G C 1 (f2 x [0,T)) satisfying 

nf' £ (-,t)=^](-,t), 

such that ?/> £ (-,0) = u (-). □ 

Theorem 5.2 (Existence and uniqueness, (3 = 0) 

and .F satisfy (1.3) and /3 = 0, traen i/iere exists a unique bounded continuous solution u°' £ on 
Tl £ x [0,T) o/tfte problem (1.4)-(1.5). 

Proof. We consider the solution u^ ,£ given by Theorem 5.1 for the choice no = Uq. Let 

u = limsup *n /3 ' e and u = liminf ^u^' 6 . 
0^0 P->o 
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Then, using Propositon 4.5, we obtain 



u < u < u < u + and \u\, |u| < HwoIIl 00 + ^ll-^lli 00 - (5-1) 

Using standard arguments similar to those in Step 1.1 of the proof of Theorem 1.1 in Section 6, we 
can show that u (resp. u) is a super- (resp. sub-) solution of the problem (1.4). The only difficulty is 
to recover the viscosity inequality on Q £ x {0}. This last inequality follows from the fact that u^' £ is 
a classical solution and then satisfies the equation also at t = 0. Finally, using (5.1), we deduce that 
u (resp. it) is a super- (resp. sub-) solution of the problem (1.4)-(1.5). Then the comparison principle 
(Theorem 3.4) implies that u < u. So u = u =: u°' £ is a continuous bounded solution of (1.4)-(1.5). 
The uniqueness of this solution follows again from the comparison principle. □ 

As a conclusion, there exists a unique (viscosity) solution u £ G C°(£f x [0, T)) of problem (1.4)- 
(1.5) defined by: 

' i/' e if > 0, 

(5.2) 

-u°' e if /3 = 0. 



u £ 



Remark 5.3 The existence of a solution could also be proven by Perron's method. 

6 Proof of Theorem 1.1 

This section is devoted to the proof of Theorem 1.1. Note that in the sequel, we use the following 
notation: 

Q T := n x (0,T), d l n T := dn x (0,T), H r :=0 x [0,T), 

and 

n £ T -.= n £ x (o,t), d l n £ T ■- dn £ x (o,r), Tf T -.= Tl £ x [o,t), 

where d l denotes the lateral boundary. 

Proof of Theorem 1.1. The proof is divided into two steps. 

Step 1: u and u are, respectively, sub- and super-solution of (1.1)-(1.2). 

Let u £ be the bounded viscosity solution of the problem (1.4)-(1.5). We define, for (x,t) £ the 
functions u and u as follows: 

u(x,t) := limsup u £ (y,s) and u(x,t) := liminf u £ (y,s). (6.1) 

y ^ x,s ^t,e ^ y^x 7 s^t,e^0 

yeTT,s€[0,T) yen e ,se[0,T) 
We start by showing that u is a viscosity sub-solution of (1.1)-(1.2). 
Step 1.1: Proof of 1 (ii) in Definition 3.5. 

Let ip € C 2 (Qt) and Pq = (a^o^o) £ such that u — (p has a zero local maximum at Pq. Without 
loss of generality, we can assume that this maximum is global and strict. Therefore we have 

Vr > 0, 35 = 5{r) > such that ip-u>6>0 in U T \B r (P ), (6.2) 
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where P r (Po) C R n x R is the open ball of radius r and of center Pq. Now, let w £ := ip — u e for some 
e > 0. Then, using the definition of u, and inequality (6.2), we infer that 

w E >- in n £ T \B r (P ), (6.3) 

for all r > and for e > small enough. Using [4, Lemma 4.2], it is then classical to see that there 
exists a sequence P* = (x*,t*) G $f T n B r (Po) such that, as e — > 0, we have 

P* — > Pq, u £ (P*) — > u(Pq) and u £ — (p has a local maximum at P*. (6-4) 

Two cases are then considered. 

Case 1 : Pq G f2 x [0,T). We argue by contradiction. Assume that there exists a positive constant 
7 > such that 



P(pt(P ) = Avp(P ) + 7 



for P G x (0, T) if /3 > 

for P 6Ox[0,T) if /3 = 1 j 



Moreover, using (6.4) and the fact that u £ is a viscosity sub-solution of problem (1.4)-(1.5), we 
conclude that 

M(P £ *)<A e M(P £ *), (6.6) 
and by Taylor's expansion, this inequality (6.6) implies 

fo t (P e *)<A V (P E *)+o £ (l). (6.7) 

Combining (6.5) and (6.7) yields: 

7 < - MPo)) + Av9(P e *) - Aip(P ) + o £ (l), 

where the right hand side goes to zero as e — > 0. This contradicts the fact that 7 > 0. 

Case 2 : Pq G 8 1 Qt- We repeat similar arguments as in Case 1. Suppose that 

min IptptiPo) - Aip(P ) , p t (P ) - F(<p(P Q )) - J^o)} = 7i > 0. (6.8) 
Inequality (6.8) implies 

P<Pt(Pb)- Av?(Po) >7i, (6-9) 

and 

^(P ) - P(^(P )) - ^(^0) > 71- (6-10) 

On the one hand, if P* G Qfn, then, using (6.7) and (6.9), we obtain a contradiction by using the 
same reasoning as in Case 1. On the other hand, if P* G d l Qj,, then using (6.4) and the fact that u e 
is a subsolution of (1.4)-(1.5) we obtain: 

^(p;)<pMp;)) + p £ M(p;), 

and then, using Taylor's expansion, we obtain 
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Finally, subtracting this inequality from (6.10), we conclude, after passing to the limit as e — > 0, that 
7i < 0; contradiction. 



Step 1.2 : Proof of l(i) in Definition 3.5. 

From Propositions 4.5 and 4.6, we can pass to the limit and get from the barriers 



tc p 

C{t + ^l + x n -l) 



if /3>0 

if /5 = and uq 



\u(x, t) — uq(x)\ < 
Therefore for t = we recover l(i) in Definition 3.5 for all (3 > 0. 



for all (x, t) G Qt- 



Step 2: Existence and convergence. 

From Step 1 we conclude that u is a viscosity sub-solution of (1. 1)-(1 .2) and by a similar manner, we 
can show that u is a viscosity super-solution of (1.1)-(1.2). Moreover u and u are bounded, because 
Proposition 4.5 implies 

I^Ij ImI ^ ll n ol|L°° + £||P||z,°° 

Then by the comparison principle for problem (1.1)-(1.2) (Theorem 3.6), we have u < u. On the 
other hand, by the definition of tt and u, we have u < u. As a consequence, we deduce, for (x, t) G ftx, 
that: 

u(x, t) = u(x, t) = lim u £ (y, s) =: u°(x,t) (6-H) 

e— >0,y— >x,s— >t 

where u° is then a continuous viscosity solution of problem (1.1)-(1.2). Moreover u° is unique, still by 
the comparison principle. Finally, we also observe that the convergence u e — > u° as e — > is locally 
uniform in the following sense: 

\\u £ -u°\\ Loa{KnWT) — ► as e^O, 

for any compact set K C £It- Indeed, suppose that there exists 9 > such that for all k > there 
exists £fc satisfying < < ^ with 

||u £fe - W°||ioc.(^ n (Q e fe x [o,T))) > ^ 1 

Then, there exists a sequence Pk £ K n (fl £k x [0,T)) such that 

|n £ (P fc )- U °(P fc )l >«• (6-12) 

Now, since KC\($f k x [0, T)) C K with compact, there exists a subsequence, denoted for simplicity 
by Pfc, such that Pfc — >■ Poo as k — > oo. Finally, taking the liminffc_ i . 00] p fe -^p 00 in the inequality (6.12) 
and using (6.11), we obtain > 9 which gives a contradiction. □ 



6.1 Other barriers and comments on another possible approach for (3 = 

Let us notice that for (3 = 0, we have natural barriers for the e-problem, which are: 

Uq' £ (x) ± C £ t 

where Uq' £ = Uq is the discrete harmonic extension (we make its dependence explicit on e) associated 
to the operator A e , and with the constant C £ satisfying 

C e >DF[v!S*\ + \\F\\ loo . 
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Then another possible approach to show the convergence of u £ to u°, could be to control the solution 
as e goes to zero, showing that: 

1) the constant C e can be taken independent on e (using the regularity of uq on dQ). 

2) the discrete harmonic extension Uq ,£ converges to the continuous harmonic extension Uq as e tends 
to zero. 

Then we could also introduce a (more classical) notion of viscosity solution in the case (3 = 0, assuming 
that at t = 0, we can compare the sub/super-solution to the initial data (taken to be equal to the 
harmonic extension Uq ,£ for the e-problem and Uq for the limit problem). Nevertheless, this other 
approach would require some additional work (to show 1) and 2)), and would not simplify the proofs. 

Remark 6.1 (Convergence of the discrete harmonic extension to the continuous har- 
monic extension) 

Notice that point 2) is a consequence of our convergence theorem (Theorem 1.1) in the case (3 = 0. 
Indeed, in the case (3 = 0, the bounded solutions satisfy 



and then 

On the other hand, the functions 



\u £ (x,t) - U q' £ (x)\ < C £ t 



u £ (x, 0) = Uq ,£ (x). 



u c (x) ± C t 

are barriers for the limit problem if 

Co > 11^ — + ll-^IU 00 - 

on 

Then bounded solutions of the continuous problem satisfy 

\u°(x,t) -«§(x)| < C t 

and then 

u°(x,0) = ug(x). 

Finally from the (locally uniform) convergence of u £ to u° in particular for t = 0, we deduce that u^' 6 
converges, locally uniformly, to Uq. 

7 Proof of Theorem 3.4 

In order to emphasize the main points, we perform the proof in several steps. 

Step 1: Rescaling. We want to reduce the problem (1.4)-(1.5) to the case e = 1 and with a 
nonlinearity F replaced by a monotone one. To this end, we introduce the new functions 

u(x,t) := e~ xt u(ex,et), v(x, t) := e~ xt v(ex, et), (x, t) G x [0,T), 

where T = T/e and A > is a constant to be determined later. We see easily that u (resp. v) is a 
sub-solution (resp. super-solution) of the following problem 

Put = A^u] - ~p\u in Q} x (0,T), 

(7.1) 

u t = F(u, t) + D 1 [u] on dQ 1 x (0, T) , 
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and 



u(x, 0) = uo(ex) for 

with 



x e tt if /3 > 

iGSfi 1 if /3 = 0, 



/3 = ep and F(u, t) = ee~ xt F(e xt u) - Xu. 
We argue by contradiction assuming that 

M := sup (u — v) > 0. 
From the definition of M, there exists a sequence 

p fc = ( x k ,t k ) eft x [0, T) such that u(P k 
M > as k — > oo. Let us choose an index ko such that 

u(P fc °) -v{P k °) > y. 

At this stage, if we write x = (x' , x n ) with x' = (x±, . . . , x n _i), we define 

I 1 2 

:= [t ~ S| + + a\x'\ 2 + 7 V1 + x n , (x,t,s) eSi'x [0,T) 2 , 

2o T — t 

where 5, rj, a, 7 > will be chosen later, and we consider 

M : = M,5 iJ))Qj7 := sup («(x, i) — u(cc, s) — ^(x, i, s)). 

Step 2: A priori estimates. Now, if we choose rj, a, 7 such that: 



T _ /A;o +«i(^°)r+7V i +^< T) 

we conclude that 

M > u(P ko ) - v(P ko ) — fy(x k ° ,t k ° ,t k °) > — > 0. 

Moreover, from the definition of M, there exists (x, t, s) G f^ 1 x [0, T) 2 such that 

«(x, f) - U(x, s) - *(x, t, s) = M > 

and 

l*-«| 2 , V , _l-/|2 



+ + a|x'| 2 + 7VI + < C 



25 T-t 
where 

C = IKHioo^xtQ,?)) + H^IL°°(n 1 x[0,T))' 

Step 3: Getting contradiction. In this step, we have to distinguish two cases: 

Case 1. (i > and s > 0) 

In this case, we have two sub-cases: 

(i) If x 6 O 1 . From the definition of (x, i, s), we see that 

u(x,t) < <^{x,t) := ~M + v(x,s) + V(x,t,s). 
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Even if u has not the required regularity, we see by a simple approximation argument that we can 
choose 

ip u (x, t) = u(x, t) for x/f, 
Then ip u is a test function for u at (x, t) and we deduce the following viscosity inequality 



_ ? ? , / Air 

(t - ty 



/3^— +/3^— ^ < A 1 [u] (x, t) - (3\u(x, t) . (7.6) 



Similarly, we have 

v(x,s) > (p"(x,s) := -M + u(x,t) - ^(x,t,s), 

and we choose 

ip v (x,s) =v(x,s) for x^x. 
We then get the viscosity inequality 

~/3^—-^- > A l [v] (x, s) - ~/3\v(x, s) . (7.7) 
o 

Subtracting the two viscosity inequalities and setting w{x) = u(x,t) —v(x,s), we get 



0<(3^— < A'lw^x) - f3\w(x) 



_']_ ^ Air 
T-tf 

< A 1 [M + ^(-,t,s)](x) 

< a (l2/'| 2 -|^| 2 ) + 7 Yl ( V 1 + Vn ~ VTT^) =:A + B, 

\y—x\=l,y£Z n \y—x\=l,y 



where in the second line we have used the fact that w{y) — ^(y, t, s) < M for all y G Q 1 with equality 
for y = x. Remark that, B = "fC(%n) with, for a > 1, 



C(o) = /(a + l) + /(a-l)-2/(a) with f(a) = VT+^, 
hence we get as in (4.12) or (4.15) 



C(a) = V t dt f dsf"{a±s)<f'{ 
_i_ jo J o 



± 

where we have used the concavity of /". Using moreover the a priori estimate (7.5) we have yj\ + x n < 
C/j, and we deduce that 

B < -C 7 4 with C' = -^-r. 
~ AC 3 

On the other hand, we have 

A = a Y {y' -x'){y' + x') < 2na(l + 2\x\) < 2n(a + 2y/aVC) , (7.8) 

\v-x\=l,i 



where we have used the a priori estimate (7.5) (q|x'| 2 < C). Finally, if we take a, 7 small enough so 
that 2n{a + 2yfa\/~C) < C'7 4 , we conclude that < A + B < 0, and hence a contradiction. 
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(ii) If x € dfl 1 . This is a similar case of the latter one where we use the same arguments. After the 
same choice of the test function, we arrive at 

+ t -^<F(u(x,t),t) + D 1 [u](x,t), (7.9) 



(t - ty s 

and 

^ > F{v(x, s),s) + D 1 [v] (x,s). (7.10) 
o 

Subtracting (7.9) and (7.10), we infer that 

< 7 (V2-i)+ E a (\y'\ 2 -\^f) + m^rt)rt)-mx,s),s)] 

\y-x\=l,y£Z n ,y n >x n 

< 7+ E a(|yf ) + [F(n(x,t)J)-F(n(xJ),s)] 

|j/-Sc|=l,I/6Z» 

+ [F(u(z,t),a) -F(w(x,s),s)]. (7.11) 

Now, if we take A > eUF'H^oo^, we get f' u < and 

F(u(x, t),s)- F(v(x, s) , s) < 0, (7.12) 

thanks to the fact that u(x,t) —v(x,s) > 0. Moreover, as \F t \ < Cq = Cq(e, A, ||i 7, ||v^ 1 .°°(iR))) we get 

\F(u(x,t),t) -F(u(x,t),s)\ < C \t-s\ < C0V2C6, (7.13) 

where we have used the a priori estimate (7.5) (\i — s\ 2 / (25) < C). Now, substituting (7.8), (7.12) 
and (7.13) into (7.11), we infer that 

< — V _. „ < 2n(a + 2^W(7) + 7 + Cov 7 ^. 
Finally, for rj > fixed, we get a contradiction by choosing 5, a and 7 small enough. 
Case 2. (i = 0ors = 0) 

For 77, a, 7 fixed, we assume that there exists a sequence 5 — > and (x, i, s) = (x 5 , F 5 , s 5 ) GO x [0, T) 2 
such that t S = or = 0. We deal with the case t S = as the other case s* 5 = is similar. Then 

«(x*,0) -77(^,5*) =M + ^(x 5 ,0,s) > — > 0. 

From (7.5), we deduce that, up to an extraction of a subsequence, we have (x s ,s s ) — > (x, 0) G 
O 1 x [0, T) as (5 — >■ 0. Therefore for /3 > we have 

< — < limsup (u(x 5 ,0) - ^(x 5 , /)) < 7i(f,0) -u(f,0) < 0, 

where we have used the comparison to the initial condition. This leads to a contradiction in the case 
(3 > or (3 = and f G dO 1 . For the case (3 = and x G O , we get a contradiction exactly as in 
case 1 i). 

It is worth noticing that, in the whole proof, we first fix r/, and then we choose respectively 7, a 
and 5 small enough. □ 
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8 Proof of Theorem 3.6 



In this section we first present the construction of an auxiliary function £ which plays a crucial role 
in the proof of Theorem 3.6. 

Lemma 8.1 (Auxiliary function) 

Let u G USC(H x [0,T)), v € LSC(Tl x [0,T)) such that u < C , v > -C , C > 0, and let 
M := sup^ x [ T Ju — v) > 0. T/ien /or a// c > 0, i/iere exists a C°° -function £ c : R n x [0, T) — >• R imi/j 
|£d < 3Co and a positive constant a c > sitc/i i/iat /or every x,y E f2 and s,t 6 [0, T) ; 

u(y, s) < £ c f-y^> ~4r~) - u ^ x ' *)' 
i/ u(x, t) — v(y, s) > M — a c , |x — y| + |i — s\ < a c and \x\,\y\ < c. 

Proof. We essentially revisit the proof of Lemma 5.2 in Barles [5]. Without loss of generality, we can 
extend the proof of our result on [0, T] by defining the end point of u and v as follows: 

u(x,T) = limsup u(y,t), v(x,T) = liminf v(y,t). 

t<T, (y,t)->(x,T) t<T, (y,i)-> (x,T) 

Let M := sup^ x j T i (u — v) with £>2c := E O; |x| < 2c} Let F = {(x,t) G l?2c x [0, T]; u(x,i) — 
i) = M}. We have u < Co, f > —Co and u - d = M on J. This implies that 

\u\,\v\ < C + M < 3C on J 7 . (8.1) 

Moreover it is easy to check that J 7 is a closed subset of £?2 C X [0, T] and that the restriction of u and 
v to J 7 are continuous. Therefore, the restriction to F of the function (x, t) i->- (ti + f)/2 is also a 
continuous function on J 7 which satisfies (8.1). We may extend this function as a continuous function 
in R n x R (still bounded by 3Co) and then, by standard regularization arguments, there exists a 
C°°-function £ such that 



~ u(x,t) +v(x,t) 
COM) 



M 

< — on J 7 and |f | < 3C on R n x R. (8.2) 
8 

In order to show the lemma, we now argue by contradiction assuming that there exist two sequences 
(xa c ,t a J, (y ac ,s a J such that for a c small enough, u(x ac ,t a J -v(y ac ,s a J > M — a c and \x ac -y a J + 

\t ac —s ac \ < cl c and such that, say, u(x ac ,t ac ) — £ ^~ ^ ac , ^ ac Sac ^ < 0. Extracting, if necessary, 

subsequences, we may assume without loss of generality that (x ac , t ac ) , (y ac , s ac ) — » (x,t). Then it 
is easy to show the convergence of u{x ac ,t ac ) — v(y ac ,s ac ) to M = u(x,t) — v(x,t). By considering 
the upper semi-continuity of u and the lower semi-continuity of v, this implies, on the one hand, 
u(x ac ,t ac ) — > u(x,t), v(y ac , s ac ) — > v(x,t). On the other hand, by using the continuity of £, we obtain 

U (x,t)-mt)= Hm («k,y-e(^,^))<o. (8.3) 

But since u(x,t) — v(x,t) = M, with (x,t) G J 7 , we deduce from (8.2) that 

~ _ s u(x,t) +v(x,t) M . _ N 3M . _ 
£(x, t) < - + — = u{x, t)- — < u(x, t), 
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which contradicts (8.3). Finally, we arrive to the result by taking £ c as the restriction of £ on [0, T) 
multiplied by a cut-off function ip £ C°°(M n ) such that ip = 1 on B c and zero outside Bi c - □ 

Proof of Theorem 3.6. The proof is divided into three steps. 

Step 1: Test function. In order to replace the nonlinearity F by a monotone one, we define the new 
functions 

u := e~ xt u and v := e~ xt v, 

where A > is a constant which will be determined later. Obviously, u (resp. v) is a sub-solution 
(resp. super-solution) of the following problem 

put = A[u] - /3\u in O x (0, T), 

du - (8 ' 4) 
u t = t— + F{u, t) on anx(0,T), 



and 



, s . . „ \ xefi if p > 

0) = u x) for ^ (8.5 

v ; v ' I x g on if /? = o, v ; 



with 

F(u, t) = e- xt F(e xt u) - An. 

Let us assume that M := sup^ x j T \(u — v) > and let us exhibit a contradiction. By the definition 
of the supremum, there exists P k ° = (x k °,t k °) 6 Q x [0, T), for some index ko, such that 

u{P k °) -v{P k °) > y. (8.6) 
Let us introduce the following constant: 

C* = IMIl°o(T2x[0,T)) + IMIl°°(?2x[0,T))- 

We now take the following notation: x = (x', x n ) with x' = (xi, . . . , x n _i), and for a, 7, 77 > (to be 
fixed later), we can approximate the functions u, v by the functions: 

£t(x, t) := u(x, i) - a\x'\ 2 - jy/l + x„ - ^ - < C*, (8.7) 



v(y, s) := v(y, s) + a\y'\ 2 + 7^ + ^ + YZ~ S - ~ C * 

and define 



M := sup (ti(x, t) — ?)(x, i)). 

xen,tG[0,T) 

In order to dedouble the variables in space and time, following the proof of Theorem 3.2 in [5], we 
define the function <6 : (U x [0, T)) 2 -)• R by: 

(x-y) 2 B(x n -y n ) 2 (t - s) 2 2(t - s)(x n - y n ) 



<£{x,t,y,s) : = ^ 1 ^ 1 ^ h 



e 2 e 2 <5 2 <5 2 

x + y t + s\ t + 
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where the parameters B,e,5 > will be chosen later. Moreover c > is a constant that will be 
defined later (only depending on a, 7, A,C*, H-Fl^oon^) and the smooth function 



£ = tc (8.9) 

is the auxiliary function associated to u and v, and given by Lemma 8.1, which shows in particular 
the following estimate 

If I < 3C* (8.10) 
We note that <j?(x, t, x, t) = 0. Then we set 

M = M B ,e,s-= _sup (u(x, t) — v(y, s) — t, y, s)) > M. (8.11) 

x,y<=fl, t,s£[0,T) 

Step 2: A priori estimates. By choosing 77, a, 7 small enough such that 



T -t k ° 



+ a|(x fc0 )r + 7V 1 + ^° < ~o~' 



it follows from (8.6) that 



M>M> u(P ko ) - v{P ko ) > ^ > 0. (8.12) 



Hence, from the definition of M, there exist sequences x k ,y k € f2, t k ,s k £ [0, T) such that 

u(x k ,t k ) - v(y k ,s k ) - &(x k ,t k ,y k ,s k ) — > M > as k -> 00, 
and, by taking A; large enough, we deduce that 



n(x fc ,t fc ) - <)(y fc , s fc ) - l»(x fc ,t,y k ,s k ) > y > 0. (8.13) 



From (8.10), it follows that 

_ f k + Ti k t k 4- s fc 7 fc + "s fc 

TO^P", — — ^)l <Ci = CtfA, C„ ||F|| LO 
Then if we take e < 5 < 1, Young's inequality yields 

T k +7i k 7 fc + s fc t k + -? k (T k -7i k ) 2 r 2 F 2 (r k -n k ) 2 r 2 s 2 

1(3* - ^)^e(^^, 4^), ^)l < ( "" £ / J + < K / J + (8.14) 

and 

|2(i* - S fc )(*£ - tn)\ ^ Ht k ~ S k ? , „ (S* ~ ^) 2 

P — + 2 — 7 2 — • (8 ' 15) 

Using (8.14), (8.15), we deduce that for B > 3 

- fk _ fc _ fc , > (x - y? (B - 3)(x B - ^) 2 (t - s) 2 C 2 <5 2 C 2 <5 2 
lX ' ,V ,S ) ~ e 2 + e 2 + 25 2 4 ~ 4 ' 

Using (8.13), we conclude that 

+ + a(|W + |(y*)'| 2 ) + + + ^ + ^ < a„ 
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with (for S < 1) 

C** = C* H — — = C#*(A, C*, ||i^||^oo (r)). 
Then, up to an extraction of a subsequence, we have 

{x k ,t k ,y k ,s k ) -> (x,t,y,s) £ (H) 2 x [0,T) 2 as fc -> oo, 

with 

_ (j 2 5 2 

u(x,t) - v(y,s) - $(x,t,y,s) = M > and $(x,t,y,s) > ^— (8.16) 

and the same a priori estimate 

+ + <*(\x'\ 2 + |yf) + 7 ^ + 7 + ^Zj + j^r= < C**. (8.17) 

Step 3: Getting contradiction. In order to get a contradiction, we have to distinguish two cases: 
Case 1. (i > and s > 0) 

(i) If x £ 00, or y £ d£l. We only study the case x £ d£l as it is similar for y £ d£l. Let 
/iflxfO.Tj-S'lbe the function defined by 

/|2 , . / -, , - , V 



(f u (x,t) :=M + v{y,s) + <$>(x,t,y,s) + a|x'| 2 + jy/1 + x n + 



T-t 

Using (8.16), we see that ip u is a test function for II and hence we obtain 

min |m(P ) - A^(Pq) - /3AH , #(P„) - W(Po),t) - |^(^o)| < 0, 

where Po = {x,t). We note that the case /3<^™(Po) — A99 u (Po) — /?AH < can be treated as in the 
sub-case (ii) below, so it is sufficient to study the case where ipf(Po) — F(ip u (Po),t) — f^-(Po) < 0, 
which implies, with x n = 0: 

11 < 2( " J 7 +F^f^,!±iy*±i > \-w.M 



(T-i)2 £2 2^1+^ 

2(1 + P)(^ W (gn-yj ^ t' i (^"-&J rEV / (g» ~ grJ fry ^ n ffi 1 ^ 
^2 2 2 v^J«-?z 2^ INt < U, (8.18) 

where 



2 ' 2 7 ' * * V V 2 ' 2 / ' 2 



and 



are the first partial derivatives of £ and P. On the one hand, from (8.16) and (8.11), we deduce that 

u(x,t)-v(y,s)>M-^. (8.21) 
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Prom (8.17) we deduce that 

\x\, \y\ < C2 = Ciip-i 7, C**) and |a? — y\ < e\J C**, \t — s\ < 5y/2C**. (8.22) 
Then in (8.9), we can choose 

c = c(a,7,A,C*, ||F||loo( R )) := C 2 

and Lemma 8.1 gives the existence of a number a c > 0. Then choosing < e < 5 small enough, we 
deduce from (8.21) and (8.22) that 

u(x, t) — v(y,s) > M — a c , \x — y \ + \i — s\ < a c and |x| , \y\ < c. 

Therefore, we deduce from Lemma 8.1 that 

f(*±* ±±l) <u(x,t)<u(x,t). 
On the other hand, for the choice A > H-^'HioWR)), we have f' u < and then we can estimate 

,tj r) ~ T ^~ t)rt) >- c \t-°\- ( 8 - 23 ) 

Now, as the terms in (8.19) and (8.20) are bounded, we conclude from (8.18), (8.23) and (8.17) that 

^-l<0(^) + 0{e) + 0{5). 

Moreover, for r/ > fixed, we get the contradiction if we take e, 5, 7 small enough so that e = e(5) < 
<5 3 < 1 and 7 < f]/T 2 . We note that in the case of (x, y) E SI x dQ, the test function for v is given by 

V v (y,s) ■= -M + u(x,t) -®(x,t,y,s) -a\y'\ 2 -J\/l + y n - jP_ ■ 

(ii) If x E 0, and y E £1. We know from (8.16) that u(x, t) —v(y,s) — &(x,t,y,s) has a local maximum 
at (x,y,t,s), where 

$(x,t,y,s) := $(x,t,y,s) + a(|x'| 2 + \y'\ 2 ) + 7 (Vl + x n + y/1 + y„) + + 

Then it is natural to apply the classical Ishii's Lemma in the elliptic case with the new coordinates 
(x = (x,t),y = (y, s)) and this is what we do. Indeed, we only use a corollary of Ishii's Lemma, namely 
Corollary 9.3 which is given in the Appendix. Applying Corollary 9.3, we get for every /x > satisfying 
IjlA < I, (with A defined in (9.2) and / is the identity matrix of M 2( - n+1 )), the existence of symmetric 
n x n matrices X, Y such that 

(ri.pi, * ) ) E V + u(x,t), (t 2 ,P2, * ) ) G V~v(y,s), (8.24) 

and 

-(£ + \\A\\)I < ( * _° y ) < A + 2/x(A 2 + A x ■ AT), (8.25) 
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where I is the identity matrix of R 2n and where all the quantities are computed at the point (x, t, y, s): 
n ■= d t $, pi := D x $, t 2 := -<9 S $, p 2 := -D y $, and 

, D 2 XX ® D 2 § \ __ ( D 2 xt <f> D 2 XS $ 

{ D 2 yx * ) and M ~ \ 



A simple computation shows that 
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Here / is the identity matrix of 1" and for all i,j G {1, n}, (I n )ij = 1 if « = j = n and otherwise, 
and (I')i,j '■= 1 if i = j < n — 1 and otherwise. Moreover 

F(z,y):=F(£(^,^),^), P' := D X ,F = D y ,F , P n := D Xn F = D y J , 

for all x = (x',x n ) S R 71-1 x R and y = (y',y n ) E R" --1 x R. Then we can write the viscosity 
inequalities for the limit sub/superdifferentials T> v(jj,s) and T> u(x,t). This gives 



and 



Pn + f3Xu(x, t) - tr(X) < 
/3r 2 + /3AU(y,s) -tr(y) > 0, 



where tr is the trace of the appropriate matrix. Substracting these viscosity inequalities and using 
(8.16), we get 

P{dt® + d s $) + /3A($ + M) < tr(X - Y), 



i.e. 



V 1} 



(T-t) 2 ' {T-s) 2 
Moreover, using the definition of 3>, we conclude that 



(x n - y n )({F)'t + {F)uQ + PK* + M) < tr{X - Y). 



P 



-{^ZW + (T Z =)2 ~ fen ~ ^n)((i ? ); + (^)UO ) + >A(-V- ("„ -77„)F) 



+ /3A 



~ l/n) 2 g ~ g) 2 2^ - S)(x n - y n ) 



<5 2 



<5 2 



< tr(X - y). 

Using the fact that S > 1 and M > 0, we deduce that 



(8.26) 



J"2 

From (8.22), we deduce that 



2/3^ - (3(x n - y n ) {(F)' t + (F)'Jt + AF) < tr(X - F). 



2/3^<tr(X-y) + 0(e). 



(8.27) 
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On (he oilier hand. Inking the matrix inequality ($.25) on (he vectors ( f J and ( . j. we gel 



with 



and 



£ \ f /|2 , A. AT\ ( £ \ ^ ..a! ^2 , 4. . aT\\ \t\2 



Using again the a priori estimate (8.17), we have 



7f 1 + ,tt4^<-CV, 



4 I (1 +3;„) 3 / 2 (l + y n ) 3 ' 2 
and we conclude that 

tr(X -Y)< 4a(n - 1) - C'7 4 - 4{x n - y n )tr(D 2 F) + 4fin\\A 2 + ■ A[ ||. (8.28) 
From (8.27) and the bounds on F in W 2,OD (W), we deduce that 

2/3^ < 4 «(™ - 1) - C7 4 + 0(e) + 4/in||i 2 + A x • AT\\. 
T 2 

Because all the terms are independent on /i, we can first take the limit fj, — > 0. Recall that (3 > 0, 
and then we can not use rj to get a contradiction as usual. We then have to get a contradiction only 
using the following inequality: 

0<4a(n-l)-C 7 4 + 0(e). 

To this end, for given 7 > 0, we get a contradiction choosing a, e small enough such that e <C 7 and 
8a(n - 1) < C'7 4 . 

Case 2. (i = 0ors = 0) 

By fixing the parameters 7, r/ and a, we assume that there exists a sequence <5 — > and (x,t,y,s) = 
(x 5 ,t ,y s ,s 5 ) G (f2 x [0, T)) 2 such that s* 5 = or i = 0. It is sufficient to study the case t = as the 
other case can be deduced similarly. Thus 

u(x s ,0) -v(y 5 ,s 5 ) = M + ${x s ,0,y 5 ,s 5 ), 

and we deduce from (8.17) (with e < <5 3 < 1) that, up to extraction of a subsequence, we have 
(x 5 ,0,y s ,s s ) -> (1,0,37,0) € (IT x [0,T)) 2 as (5 0. Hence in the case /3 > 0, we conclude that (using 
(8.16)) 

< M < limsup (M + $(x s ,0,y 5 ,s 5 )) = limsup (u(x 5 ,0) - v{y s , s 5 )) < u(f,0) - v(W,0) < 0, 
<5->0 5^0 

where we have used the comparison to the initial condition. This gives a contradiction in the case 
/3 > and in the case /3 = if x G <90 . In the case /? = and x £ £1, we get a contradiction exactly 
as in the case 1 i). 

In the whole proof, we first fix r\ and B. After that, we fix respectively 7, a, 5, e (and fi). □ 
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9 Appendix 



This appendix is dedicated to the proof of a corollary of Ishii's lemma (Corollary 9.3), which is used 
in the proof of Theorem 3.6. First, we recall the elliptic sub and superdifferentials of semi-continuous 
functions and the classical Ishii's Lemma. In all what follows, we denote by S n the set of symmetric 
n x n matrices. 

Definition 9.1 (Elliptic sub and superdifferential of order two) 

Let U be a locally compact subset o/R n and u G USC(U). Then the superdifferential T> + u of order 
two of the function u is defined by: (p,X) G R n x § n belongs to T> + u(x) if x G U and 

u(y) < u(x) + (p, y - x) + ^(X(y - x), y - x) + o(\y - x\ 2 ) 

as U 3 y — >■ x. In a similar way, we define the subdijferential of order two by T>~u = —T> + (—u). We 
also define: 

( (p, X) G M n x S n , 3(x n , Pn , X n ) G U x R n x S" 
V u(x) = < such that (p n ,X n ) G V + u(x n ) 

[ and (x n ,u(x n ),p n ,X n ) ->■ (x,u(x),p,X) 

The set D u(x) is defined in a similar way. 

Now recall the classical elliptic version of Ishii's Lemma. 

Lemma 9.2 (Elliptic version of Ishii's lemma) 

Let U and V be a locally compact subsets ofW 1 , u G USC(U) and v G LSC(V). Let ip : U x V ^ R 
be of class C 2 . Assume that (x, y) i-> u{x) — v{y) — ip(x, y) reaches a local maximum at (x,y) G U x V. 
We note p\ = D x (p(x,y), p2 = —D y (p(x,y) and A = D 2 (p(x,y). Then, for every fi > such that 
fj,A < /, there exists X, Y G §> n such that: 



( Pl ,X)eV^u(x), (p 2 ,Y)eV v(y) 

,1 
/' 



(- + llW* _°y)<A + ^, < M > 



where I is the identity matrix ofR 2n . The norm of the symmetric matrix A used in (9.1) is 

\\A\\ = sup{|A| : X is an eigenvalue of A} = sup{| < A£,£ > | : |£| < 1}. 

For the proof, we refer the reader to Theorem 3.2 in the User's Guide [15]. 
Then we have: 

Corollary 9.3 (Consequence of Ishii's lemma) 

Given T > 0. Let U and V be a locally compact subsets of R n , u G USC(U x [0, T)) and v G 
LSC(V x [0,T)). Let $ : U x [0,T) x V x [0,T) -)• R be of class C 2 . Assume that (x,t,y,s) h-> 
u(x,t) — v(y,s) — &(x,t,y,s) reaches a local maximum in (x,t,y,s) G Ux [0, T) x V x [0, T). Computing 
the following quantities at the point (x,t,y,s), we set t\ = dt&, T2 = —d s Q, p\ = D x &, p2 = —D y Q 
and 

1 Aj A 2 J V D yx® D 2 yy § )> V D yt® D%s® J \ D 2 st $ D 2 S $ 

(9.2) 
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Let I be the identity matrix o/M 2 ( n+1 ) . Then for every fi > such that \xA < I, there exists X,Y G S r 
such that (where we denote by * some elements that we do not precise): 



X * \ + - / Y * 

(n,Pi,( + ^ J) u{x,t), (t 2 ,p2, ( ^ + 



G X> w(y,s), 



and 



+ < ( * _° y j < A + 2^ + A! • A?'). 



where I is the identity matrix of. 



p2n 



Proof. Because $ £ C 2 and u(x,t) — v(y,s) — &(x,t,y, s) admits a local maximum in (x,t,y,s), 
we can then apply the elliptic Ishii's Lemma (Lemma 9.2) with the new variables x := (x,t) and 
y := (y,s). We obtain, for every /i satisfying pA < I with the 2(n + 1) x 2(n + 1) matrix A = D 2 $, 
that there exists IT 6 S n+1 with 



X = 
such that: 

and 



X C 
l C 5 



Y 



Y D 



with X,Y£§ n , C,D£R n and 5,p£ 



l D p 

(pi,ti),x) £v + u {x,t), ((p 2 ,r 2 ),y) er«(j,s) 



-(^+IIAII)/< 



V 



*c 


c 

5 








-y 


-p , 



\ 



< A + /^A 2 



(9.3) 



We first remark that the matrix A (defined in (9.2)) is obtained from A by relabeling the vectors of 
the basis (going from coordinates (x, t, y, s) for A to coordinates (x, y, t, s) for A). Therefore we have 
\\A\\ = \\A\\ and the condition pA < I is equivalent to pA < I. 

Next, for G M ra , applying the vector V = 0, n, 0) T to the matrix inequality (9.3), yields with 
V = (£,r ] ) T : 



-(l + \\A\\)<i.v,v> < <(q _° y )-^> 

< < (A + pA 2 ) -V,V > =:J. 
Remark that using the relabeling of the vectors of the basis, we have 



(9.4) 



< A k V, V >=< A k V, V > with V 



tt V 



for k = 1,2 



We compute 



V 




A Ai 
Af A 2 







y T Ay 



and 



V 




A Ai 
Af A 2 



W = y T A 2 y + y T AiAf v + 2 < Ay, a\v >< 2 (|Ay| 2 + |Af v\ 2 ) 
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This gives 



J < V T AV + 2fi (\AV\ 2 + \AfV\ 2 



) 



and then (9.4) implies 



(7 + < 



( 



X 
-Y 



) 



< A + 2fi(A 2 + A x ■ A\) 



which achieves the proof. 



□ 
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